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Abstract 

We develop an adequate description of non-topological solitons with a small charge, for 
which the thin-wall approximation is not valid. There is no classical lower limit on the charge 
of a stable Q-ball. We examine the parameters of these small-charge solitons and discuss the 
limits of applicability of the semilcassical approximation. 
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Non-topological solitons 0, in particular Q-balls 0, are known to exist in the limit of a 



large charge. Modern theories often envision the physics beyond the Standard Model (SM) as 
associated with a number of scalar fields with various charges. All supersymmetric generaliza- 
tions of the SM, in particular the minimal model, MSSM, allow for a variety of baryonic and 
leptonic balls built of squarks and sleptons, and the Higgs scalars ||. These objects can be 
produced in the early Universe and can lead to interesting cosmological consequences. Q-balls 
with a small charge can be produced more easily at high temperatures. It is of great interest, 
therefore, to understand how small the charge of a stable Q-ball can be, and what are the 
parameters of such solitons in the limit of small Q. 

The usual description of Q-balls || relies on the so-called thin-wall approximation and 
is only valid for a very large charge Q. A naive extrapolation of these results beyond their 
domain of validity would lead one to conclude that no stable low-Q solitons can exist. We will 
prove this not to be true, at least as long as the semiclassical treatment of Q-balls remains 
appropriate. There is no classical lower limit on the charge of a stable Q-ball. We will show 
that very small Q-balls (Q-beads) with charges Q > 1 can exist and we will develop a formalism 
that yields an adequate description of these solitons. 

1 Abelian Q-balls 

Let us consider a field theory with a scalar potential U(ip) which has a global minimum 
£7(0) = at ip = 0. Let U(ip) have an unbroken global[| U(l) symmetry at the origin, ip = 0. 
And let the scalar field ip have a unit charge with respect to this U(l). 

The charge (taken to be positive for definiteness) of some field configuration <p(x,t) is 



Since a trivial configuration <p(x) = has zero charge, the solution that minimizes the 
energy 



1 Q-balls associated with a local symmetry have been constructed Q]. An important qualitative difference 
is that, in the case of a local symmetry, there is an upper limit on the charge of a stable Q ball. 
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and has a given charge Q > must differ from zero in some (finite) domain. We will use 
the method of Lagrange multipliers to look for the minimum of E at fixed Q. We want to 
minimize 



E u = E + u 



(3) 



where to is a Lagrange multiplier. Variations of <f{x,t) and those of u can now be treated 
independently, the usual advantage of the Lagrange method. 
One can re-write equation (|]) as 



= J d 3 x -\d t (p — iui(p\ 2 + J d" 
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where 



um = u{<p) 
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We are looking for a solution that extremizes S^, while all the physical quantities, including 
the energy, E, are time-independent. Only the first term in equation @ appears to depend on 
time explicitly, but it vanishes at the minimum. To minimize this contribution to the energy, 
one must choose, therefore, 

<p(x,t) = e^vix), (6) 
where <p(x) is real and independent of time. For this solution, equation (HD yields 



Q = uj J ip 2 (x) d 3 x 
It remains to find an extremum of the functional 



(7) 



d 3 x 



-\V V (x)\ 2 + U^(x)) 



+ ujQ, 



with respect to uj and the variations of <p(x) independently. We can first minimize £ u for a 
fixed uj, while varying the shape of <f(x). This, however, is identical to the problem of finding 
the bounce <fu}(x) for tunneling in d = 3 Euclidean dimensions |^, ^ in the potential U u (ip) 
(Fig. 1). The first term in equation (]8|) is then nothing but the three-dimensional Euclidean 
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Figure 1: Finding a Q-ball is equivalent to finding a bounce that describes tunneling in the 
potential U u (ip) = U(cp) — (l/2)u 2 {p 2 . Thin-wall approximation is good for large Q (thin solid 
line), but breaks down when Q is small and, therefore, to is almost as large as the mass term at 
the origin. In the latter case (dashed line), the "escape point", <p(0) is close to the zero of the 
potential and is far from the global minimum. This allows for an alternative approximation 
that is good outside the thin-wall regime. 

action S3 [(/^(rc)] of this bounce solution. In what follows we will use this analogy extensively. 
Over the years, the problem of tunneling has been studied intensely, and many properties of 
<p u (x) and Sz[{pu(x)] are well-known. 

If the following condition is satisfied, 

U ((f) J if 2 = min, for ip = ip Q > (9) 

the corresponding effective potential Uu (<p), where u = ^2U (ipD/tp 2 , will have two degen- 
erate minima, at ip = and ip = ip . 

The existence of the bounce solution (p w (x) for loq < uj < U"(0) follows 0, 0] from the fact 
that U w (<p) has a negative global minimum in addition to the local minimum at the origin. 
From Ref. [|7| we know that the solution is spherically symmetric: <p(x) = <p(r), r = 

The soliton we want to construct is precisely this bounce for the right choice of u, namely 
that which minimizes S w . The last step is to find an extremum of 

S u = S 3 [<p IJ (x)]+uQ (10) 
with respect to u. The conditions for the existence of such an extremum will be discussed 
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Figure 2: Solitons with a smaller charge Q have larger values of uj. The curved dotted 
line shows S , 3 [(^ tJ ( function of uj. The punctured straight lines 1 and 2 correspond to 

ujQi and for Q\ > Qi- The minimum of S^ip^x)] + is reached at uj = c^i, or dj 2 , 
respectively. If Qi > Q2, then Wi < uj 2 . 

below. The values of uj range from ujq, the minimal value for which U w (<p) is not everywhere 
positive and, therefore, the bounce solution can exist, to uj = U"(0), at which point Uui^) has 
no barrier. The latter does not automatically mean that a non-trivial solution of the equations 
of motion does not exist ||. However, for our purposes, we will only have to consider uj in 
the range uj < uo < U"(0). 

In the limit u — >• Uq + 0, the bounce (and, therefore, a Q-ball) solution can be analyzed in 
the thin- wall approximation 0. However, for larger values of u, the thin- wall approximation 
breaks down and so does the existence proof of Ref. 0. In our analysis, we will use a different 
approximation ||, that which is valid in the limit of large u. 

Small, near-critical values of uj ~ uj , correspond to the large values of charge Q. There 
is a simple way to understand this. The first term in equation (|10|), S^ip^x)), is a monotone 
decreasing function of uj. (The smaller the barrier, the more probable the tunneling is. Thus 
S3 must decrease with uj.) However, the last term in (|3|), ujQ, increases with uj. Therefore, if 
there is a minimum, it will be achieved for a smaller value of uj if the Q is larger, as illustrated 
in Fig. 2. 

If, however, Q is sufficiently small, then the value of uj that minimizes 8 U is large enough 
to destroy the near degeneracy of the two minima (Fig. 1), and the thin- wall approximation 
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breaks down. 



2 Thick-wall approximation 

There is, however, a powerful analytical approximation^ || that can be used for calculating 
S3 [<^(x)] in the limit of very non-degenerate minima. It is based on the observation that for 
large u (this case is analogous to tunneling into a very deep minimum), the "escape point" 
92(0) (which is also the maximal value of (p inside the Q-ball) is close to a zero of U w ((p), and 
is far from its minimum. As the barrier becomes smaller, the escape point <^(0) moves closer 
to the origin (Fig. 1). In this limit, one can neglect the dynamics at large (p and retain only 
the quadratic and cubic terms in the potential. 

We will apply this strategy to minimize of 8 U in equation (f|). Let us consider a potential^ 
U(ip) = \M 2 ip 2 — Aip 3 + A4V? 4 , where <p has a unit charge with respect to some global U(l) 
symmetry unbroken at tp = 0. Then 

u u (<p) = ^(m 2 - w V - V + W (11) 

The bounce for the potential U^ip) for large uj (that is < (M — u>)/(M + uj) <C 1) is 
the same as that for A4 — > 0. As we show below, the condition for the quartic terms to be 
negligible, A^ 4 <C Aip 3 when (1/2)(M 2 — uj 2 )ip 2 ~ Aip 3 , is satisfied whenever the charge is 
small enough. One can, therefore, neglect the quartic term in equation (|TT|) and introduce 
some new dimensionless variables for the space-time coordinates, = (M 2 — and for 

the dynamical field, ip = ipA/(M 2 — uo 2 ). In terms of these variables, the expression for S w 
becomes 



+ wQ. (12) 



The first term in equation (|T2]) is a dimensionful coefficient times the action of the 
bounce in the potential Y(ip) = (1/2)?/' — ip 3 . It has been computed || [T2|] numerically: 

2 Another approach to calculating ^3 [ip^ (x)} analytically outside the thin- wall limit was proposed in Ref. 
p0| . To calculate the bounce action numerically, one can use the Improved Action method p!l[ , which is 
particularly useful for systems with many scalar degrees of freedom. 

3 The ip 3 term represents some Z7(l)-symmetric cubic interaction, e.g., (ip^ ip)( 3 / 2 \ In the MSSM, the 
requisite cubic interactions arise from the tri- linear couplings of the Higgs field to squarks and sleptons Q . 
For clarity, we ignore the "flavor structure" of these cubic terms for now. 
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~ 4.85. Its precise value is not important, we will only need the fact that is independent 
of M, u, A, and Q. The corresponding bounce has a radius ~ 1 in dimensionless units ||. 
We have, therefore, extremized with respect to the variations of the or, equivalently, 
those of ip(x). What remains is to find a minimum of 



(M 2 - u 2 ) 3 /- 

with respect to u>, < uj < M. This is possible as long as 



£ U = S^ — VwQ (13) 



ZS^M 2 2 

Minimum is achieved at uj — uj m i n = M 
of the soliton can be expanded in powers of e: 



QA ' <i (14) 



1 1/2 

1 + v 1 — 4e 2 )/2 . The resulting emergy (mass) 



E = QM 



1 - -e 2 - -e 4 - 0(e e 
6 8 v 



(15) 



Since the mass of the soliton (|i"5|) is less than QM, it is stable with respect to decay into the 



cp quanta. The size R of the soliton is ~ 1 in dimensionless units a/ M 2 — tu 2 R, that is 



R- 



(M 2 - u 2 ) 1 / 2 » e M (l + l -e 2 + ^e 4 + 0(e 6 )) (16) 

Two conditions must be satisfied in order for the approximation we just used to be valid and 
self-consistent: (i) the quartic terms must be much smaller than the cubic and quadratic terms 
of the U in the vicinity of </?(0) (and, therefore, for every (p(x), because <p(Q) = max[ip(x)]); 
and (ii) S u in equation (|13D must have a minimum for < u < M. These two requirements 
set the following respective constraints on the charge Q. 

Q < 3£^M/VA^4~ 14.6M/VAIA 

(17) 

Q < 3S, 4} M 2 /2A 2 w 7.28 M 2 /A 2 
When these constraints are strongly violated, the thin-wall approximation can be used. 

We conclude that in the limit of small Q, the thin wall approximation must be replaced 



by the approximation (12|). The existence of the solition was proven in two steps. First, we 



noted that S w has an extremum with respect to the variations of <f(x) for fixed uj. This relies 
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on the existence proof in Ref. ||. Second, using the "thick- wall" approximation, we showed 
that £<j has a minimum for some u < M. 

The expressions for the energy fll5f) and the size (|16| ) of the soliton are accurate when the 



cubic (tri-linear) couplings in the potential are not too large, and when the charge is small 
(equatiuons QTJD and (|T7])). 

There is no classical lower limit on charge: no matter how small Q is, there is a value of 
uj close to M, for which the energy of the configuration is minimal. However, for reasons of 
quantum stability, Q must be an integer. Therefore, Q > 1. Also, as we show below, in the 
limit Q — > 1, the quantum corrections can be significant. 

3 Classical stability and quantum corrections 

We would like to examine the second variation of energy and prove the stability of the soliton 
with respect to small variations that conserve charge. In the sector of a given charge Q (the 
only subspace of the functional space in which we are interested), E and £ u coincide, and so 
do their variations. The time derivatives enter only in the first term of equation (|j), which is 
a non-negative function minimized by our solution (^]). What remains is to consider (5 2 E) Q 
or (5 2 S UJ ) Q with respect to the variations of the time-independent part of @, 6ip(x). For 
arbitrary 5<p, the variation of to is induced, so that the charge conservation constraint (0) is 
satisfied. Starting with either of the expressions (0) or (pi), one finds that 



(S 2 E) Q = J 5(p[-A + U"(<p(x)) + 3u 2 ]5 V d 3 x (18) 

where 3u 2 is the contribution of the second variation of the (l/2)u 2 (p 2 term under the constant 
charge constraint: 5 2 {(l/2)u) 2 ip 2 } Q = 3u 2 5(p 2 . If we expand Sip = J^Ci^i m terms of the 
orthonormalized eigenvectors of the differential operator in square brackets, ipi, corresponding 
to the eigenvalues Aj, then (5 2 E) Q = J^cf^i- Therefore, if the operator in square brackets in 
equation (p~8f ) has only positive eignevalues, then (5 2 E) Q is positive definite and the soliton is 
stable with respect to small perturbations. We, therefore, consider the following eigenvalue 
problem: 

[-A + f/"(v^(x)) + 3cu 2 ]^ = Xifa (19) 
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with the boundary condition ^(oo) = and the normalization condition / ipi^Pj — 

It is easy to see that for large enough u> operator ( jl9|) has positive eigenvalues only. Indeed, 
equation (|19|) is just a Schrodinger equation for the potential U"{(p u {x)) + 3u 2 . Since the 
potential U(<p) has a minimum at the origin, there exists a value <p CO n, such that U"(<p) > 
for < <p < (p con . For large u, <p(0) is small (Fig. 1) and, for a large enough to (while 
still uj < U"(0)), <p(0) < (peon- Since Vx, <p(x) < <p(0), equation ( |19|) describes a quantum- 
mechanical bound state of energy Aj in the potential that is everywhere positive. Clearly, Aj 
is then also positive (and, in fact, Aj > 3a; 2 ). 

We have proven the stability of the Q-ball in the limit of large u (small Q) with respect to 
small perturbations. Coleman proved @ that Q-balls are stable with respect to all, not only 
small, deformations in the limit of large Q. We do not have a rigorous proof of the soliton 
stability in the intermediate region, although it seems plausible. 

In order to evaluate the validity of the semiclassical approximation, one has to examine 
the magnitude of the quantum corrections to the mass of the soliton. Semiclassical results 
are reliable if the quantum fluctuations around the soliton are not large in comparison to its 
energy. The spectrum of such fluctuations is given by the same operator, 5 2 E/5(p 2 , in the 
soliton background, renormalized with respect to the oscillations around the trivial vacuum 
solution. The high-frequency modes cancel out, but the low-frequency spectrum of 5 2 E/5ip 2 
around the soliton will contain discrete levels which we would like to estimate. 

Fortunately, we know something about the low-energy modes of a different operator, which 
differs from ([19]) by a constant. Indeed, for the three-dimensional bounce solution in the 
potential fT^ (<£>), 

^ = -A + U^(x)) = -A + U u (<p u {x)) - u 2 . (20) 

The spectrum of the operator ( p0|) was studied in Ref. 0. We know that it has one 
negative eigenvalue A x < 0. However, we just proved that Ai = Ai + Au 2 > 0. Therefore, the 
lowest eigenvalue of the operator fll9|) , Ai, is in the range < Ai < Auo 2 . One can, therefore, 
expect the corrections to the soliton mass squared to be (at the most) of order uj 2 , which, in 
the limit of small Q, is ~ U"(0). These corrections will be small in comparison to the soliton 
mass squared ( |T5D if Q 2 ^> 1. However, the semiclassical approximation can become unreliable 
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for calculating the masses of the solitons with charge Q ~ 1. (At least, we don't have a proof 
to the contrary.) Nevertheless, even in this limit the size R ~ e -1 M of the soliton remains 
large in comaprison to its De Broglie wavelength, which is an indication that semiclassical 
treatment may otherwise be appropriate for Q ~ 1. 

4 Virial theorem 

We will now prove a useful virial theorem which is valid for any Q, and requires no approxi- 
mation. (In doing so, we will also illustrate that the soliton is a global minimum with respect 
to the size variations.) Let's consider a one-parameter family of functions obtained from the 
solution <p by expanding (contracting) it by a factor a: 



ip a {x) = tp(ax) (21) 



The energy of (p a is 



E ^^2]$jT X +aT + a3V - (22 > 
where T = J i(Vy?) 2 d 3 x is the gradient energy, and V = J U(ip)d 3 x is the potential energy 

of the Q-ball, both positive. 

Theorem: 

T + 3V = nr 2, 3 ( 23 ) 
2 J (p z d 6 x 

Proof. Since (p is a stationary point of E with respect to all variations, it must, in particular, 
be an extremum with respect to scaling fl2~ip . Therefore, dE a /da = at a = 1. This yields 
relation (f23|). □ 

If Q — 0, there is no non-trivial solution, since T > and V > would imply T = V = 
0. Equation (|23| ) shows, in particular, how the non-topological solitons evade the Derrick's 
theorem [|1^ . 

We note that, since the coefficients of 1/a 3 , a, and a 3 in equation ( [22] ) are all positive, 
the only exteremum of E a with respect to a (the size of the soliton) is a global minimum. 

In summary, Q-balls with a small charge exist and are classically stable. They can be 
treated semiclassically at least as long as Q 2 3> 1. For Q « 1 the quantum corrections to the 

9 



soliton mass are not necessarily small, but the Q-ball remains an extended object, whose size 
is large in comparison to its De Broglie wavelength. 

The author would like to thank M. Shaposhnikov for many interesting discussions, K. Lee 
for a number of useful comments, and L. Alvarez- Gaume, S. Dimopoulos, G. Dvali and 
G. Veneziano for helpful conversations. 
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